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Abstract 

We describe two independent methods to calculate the angular distribution of muons 
after traversing a thick scatterer due to multiple Coulomb scattering. Both methods 
take into account the nuclear size effect. We demonstrate a necessity to account for 
the nucleus extension as well as incoherent scattering on atomic electrons to describe 
the muon scattering at large angles in thick matter layers. The results of the two 
methods of calculations are in good agreement. 
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1 Introduction 



Multiple scattering of muons in the field of atoms is of interest for numer- 
ous applications related to muon transport in matter, in particular, for simula- 
tion of the response of high-energy particle detectors (including stochastic par- 
ticle deviation in magnetized steel spectrometers), radiation protection tasks 
at accelerators, muon-induced background evaluations for cosmic ray neutrino 
experiments, etc. In many cases, reliable estimations of angular and lateral 
muon distribution functions at the level of probability of about 10~ 2 — 10~ 4 
and lower are necessary. Furthermore, renewed interest to particle momentum 
evaluation on the basis of precise measurements of multiple scattering effect 
(utilising the scales and coordinate accuracies of novel muon and neutrino 
detectors) calls for adequate accuracies in the theoretical description of the 
phenomenon. 

Several multiple scattering theories have been published (see review [1]) 
which are concerned with angular distribution of particles passing without 
substantial loss of energy AE <C E using the small-angle approximation of 
scattering angle x- This approximation assumes that x is small and consists 
in replacing sin(%) tan(x) m x an d the upper limit 7r for x by infinity. 
Gaussian distribution function is widely used for various estimations, however 
it can hardly serve as a satisfactory tool for consideration of low probabil- 
ity region (large deflection angles), especially at moderate absorber thickness. 
On the other hand, Moliere theory [2], [3] of multiple scattering, the results 
of which are employed nowadays in most of the transport codes, cannot be 
used for description of high-energy muon scattering either. The reason is that, 
due to its lagre mass, muon is capable of passing large thicknesses of material, 
and hence the probability of scatters with impact parameters comparable with 
nuclear size is not small. Moliere theory (as well as other approaches devel- 
oped within the point-like nucleus approximation) heavily overestimates the 
probability of large angle deflections. 

Several calculations have been carried out to take into account the charge 
distribution in nucleus. Two methods [4] for evaluating multiple scattering on 
extended nuclei were given, one of which is of interest as a general method for 
any charge distribution. The expression for angular distribution was obtained 
in Ref. [5] under the assumption that the charge distribution in nuclei and pro- 
ton is Gaussian. Unfortunately, there are several misprints in the work [5]. 

Moliere theory regarded elastic collisions against the screened Coulomb 
field of atomic nuclei only. Collisions with Z atomic electrons also contribute 
to multiple scattering, especially in light elements. This effect is often intro- 
duced by replacing Z 2 with Z(Z + 1), which serves only to estimate the order 
of magnitude of the effect of inelastic scattering. The account of inelastic col- 
lisions in the Moliere theory has been done in Ref. [6] for electrons and heavy 
particles. 

In this paper, we present two methods for calculation of the muon angular 
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distribution function due to multiple Coulomb scattering (MCS) in matter. 
The procedures take into account the nuclear size effect and contribution of 
incoherent scattering on atomic electrons. 

The paper is organized as follows. In Sec. 2, we give the general rela- 
tions and definitions used for description of multiple scattering of muons in 
the field of atoms. A mixed Monte Carlo approach to the MCS considera- 
tion is described in Sec. 3. In Sec. 4, contribution of incoherent scattering on 
atomic electrons and the multiple scattering in mixtures and compounds are 
discussed. Basic expressions of Moliere multiple scattering theory are given in 
Sec. 5, whereas the equations of modified Moliere theory taking into account 
nuclear size effect are described in Sec. 6. At the end of this section, we com- 
pare the angular distributions of muons obtained by Monte Carlo method and 
analytically. In the conclusion, the main results of the paper are summarized. 



2 General relations and definitions 



Differential cross section (per gram of matter) for scattering in Coulomb 
field of nuclei in the small angle approximation can be written as follows [7]: 

^-*-7(3 2 f-*7(f) ! ^ W 

where N\ v is the Avogadro number, A is the atomic weight, p and v are 
momentum and velocity of the particle, respectively. In the second relation, 
classical electron radius r e and electron mass m e are utilised as natural units of 
distance and energy. Hereafter, whenever it does not cause misunderstanding, 
we assume ft — c — 1. For high-energy (E ~ p, f3 ~ 1) muon scattering on 
heavy target through small angles, the transverse momentum transferred to 
the particle q = Ex, and practically coincides with the total 3-dimensional 
momentum transfer. The cross section may be re-written as: 

da (q) = 8nN Av — {r e m e f — . (2) 



Noteworthy, in this ultrarelativistic regime the cross section depends only 
on transverse momentum and not on the particle energy. The integral cross 
section for scattering with transverse momenta greater than q is: 

a-\q) = AnN Av ^ ( r -^A 2 . (3) 
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Differential cross section for extended nuclei (taking into account also atomic 
screening) in the Born approximation may be described by means of formfac- 
tors: 

da(q) = da (q) {(F N - F a f + 1 (l - jft) F p 2 } , (4) 



where F a (q), F N (q), and F p (q) are elastic atomic, nuclear, and proton formfac- 
tors, respectively. The second term in braces represents the (order-of-magnitude) 
correction for quasi-elastic processes (with excitation or disintegration of the 
nucleus), if incoherent scattering of the projectile on individual protons inside 
nucleus occurs. For hydrogen (Z =1) this term should be omitted. 

Formfactors are normalised in such a way that F(0) = 1 and F(oo) = 0. A 
big difference between the sizes of atom and nucleus allows, as a rule, to con- 
sider the influence of atomic and nuclear formfactor separately: in the region 
where F a is essential, F^ = 1; on the contrary, when F^ appreciably deviates 
from the unit, F a is negligibly small. Hence, in a wide range of momentum 
transfers (two - three orders of magnitude, depending on the substance), and, 
as a consequence, over 4-5 decades in interaction frequency determined by 
Eq.(3), individual scatters obey a simple Rutherford law. 

The functional dependence of the formfactors is determined by the charge 
density distribution inside the object. For light and medium nuclei, a suffi- 
ciently good description may be reached on the basis of two-parameter dis- 
tribution function suggested by Fermi. However, the use of accurate charge 
distributions entails serious computational difficulties and cannot be done in 
analytical form. Therefore, reasonable approximations are often used (for ex- 
ample, Gaussian or exponential distribution of charge density). For consider- 
ation of multiple scattering process, such approximations are justified by the 
fact that at moderate transferred momenta the same expansion can be written 
for all formfactors 

F N (q)~l-(qR N ) 2 /6, (5) 



where Rn is rms radius of charge distribution. An appreciable difference be- 
tween different formfactors appears only in rare collisions with large q. Ana- 
lytically, within the framework of Born approximation [8] 



Fn(q) = exp 



{qR N y 



and F]sr(q) 



1 + 



(qR 



N 



12 



(6) 



for Gaussian and exponential charge distributions, respectively. Rms-radii for 
practically all elements (and, moreover, for numerous isotopes) can be found 
elsewhere (e.g., see the compillation in Ref.[9]). For light and medium nuclei, 
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the A-dependence of Rn may be parameterized as 

R N = 1.27A - 27 fm or q N = 1/R N = lttA' - 27 MeV. (7) 



For proton, R p = 0.85 fm, and q p = 232 MeV. 

The important parameter in consideration of multiple scattering phe- 
nomenon is a so-called characteristic value q c of transverse momentum, which 
determines on the average one scatter with q > q c at a given layer thickness 
X(g/cm 2 ): 

q 2 c =4nN Av ^-(r e m e ) 2 X. (8) 



Sometimes, the following regions are distinguished in the scattering process: 
multiple scattering (q <C q c , numerous scatters), plural scattering (several, 
but few, scatters with q ~ q c ), and single scattering (low probability scatters 
with q ^> q c ); in the latter domain, scattering angle distribution function ap- 
proaches the behaviour of the differential cross section. 

Statistical consideration of multiple scattering is appropriate when q c ^> 
1/R a , where R a is the effective radius of the atomic electron charge distribu- 
tion: 

183Z-V3 

* = 1 /*=2JIW (9) 



Hence, the applicability of theory is usually limited to thickness of layers: 
X > 4 • 1Q- 4 A/Z 4 / 3 g/cm 2 . (10) 



3 Monte Carlo approach 



The idea of MCS consideration by means of a mixed Monte Carlo technique 
(see, e.g., Ref.[10]) consists in the separation of multiple scattering through 
small angles, or with small transferred momenta (which are treated in Gaus- 
sian approximation), and of scatters in the region of large q, which are sim- 
ulated explicitely. This approach can be conveniently used for calculations 
with accurate cross section in high q region (in particular, including nuclear 
formfactors) . The boundary q^ between these two domains has to satisfy the 
following inequalities: 

Qa < Qb < Qn; Qb < q c - (11) 
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The first of these two conditions ensures Rutherford behaviour of the cross 
section around while the second one provides the regime of plural or mul- 
tiple interactions for scatters that are simulated individually. 

In order to calculate the rms-deviation for small angle scatters with q < q^, 
which are treated in a " continuous" way, the following integral has to be eval- 
uated: 

hl-FM)] 2 ^- (12) 





A similar integral comes into cross section formulae for many electromagnetic 
interaction processes (for example, radiation logarithm in the bremsstrahlung), 
and a ready result of work [11] can be used. In the frame of the Thomas- Fermi 
model, one can find for the "restricted" value of squared transverse momentum 
accumulated in the layer X of matter: 



Z 2 

(Q 2 )rstr = 8nN Av (r e m e ) 2 — X 



lnXZ-^ + mA-i 



(13) 



where i^=183. Using Eq.(8) and Eq.(9), we can re- write the latter relation as 



(Q 2 )rst r = 2g c 2 ln 



Qa 



(14) 



For light and medium elements, Thomas-Fermi model of the atom does not 
provide high accuracy. For example, for hydrogen it is necessary to take 
K =202.4; the values of K for various elements calculated with Hartree-Fock 
model may be found elsewhere [12]. 

Random scatters with q > qb are simulated by means of a usual technique, 
with the mean free path defined by the integral cross section (3). Random q q 
are sampled as 

Q, = <lb/Vs, (15) 



where q is a random number with a uniform distribution between and 1; q ( 
is accepted with a probability equal to the combined nuclear formfactor T% 
(including elastic and quasielastic terms): 

ffo) = + ! (1 - (16) 



Again, similar to Eq.(4), only the first term remains for hydrogen (with nu- 
clear formfactor equal to that of proton). 
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Although the consistency of the approach based on the partition of in- 
dividual scatters into "continuous" (q < q b ) and "discrete" (q > q b ) parts 
is almost obvious, we have performed comparative calculations for different 
thicknesses of absorber (1 cm and 100 cm iron), each with two substantially 
different values of q b (0.1 and 1 MeV for thin absorber, 1 and 10 MeV for 
thick layer); 10 6 muons for the lower value of q b and 10 8 for the higher one 
have been traced. Noteworthy, the number of simulated individually scatters 
(per muon) varies as 1/qf (see Eq.(3)), and hence it was 100 times different 
(about 15 per layer for the higher value of the threshold and ~ 1500 for the 
lower one). 

Differential distributions of events in the value of accumulated transverse 
momentum Q are shown in Figs.l and 2 for 1 cm and 100 cm iron targets, 
respectively. A big difference between calculation results for point-like and 
finite-size nucleus is obvious already for 1 cm layer. One the other hand, the 
results obtained with different boundary values q b agree within the statistical 
accuracy in the whole range of angular deviations. 

Monte Carlo approach allows to easily probe different models of the nuclear 
formfactor. In Fig.3, comparison of calculation results obtained for several ver- 
sions of the formfactor is given. Integral distributions in Q (i.e. probabilities 
of the deflection at angles greater than Q/E) after passing 100 cm iron layer 
are presented. Solid curves (bottom to top) in the figure correspond to: Gaus- 
sian charge distribution in nucleus, incoherent scattering on protons being 
neglected (i.e., the second term in Eq.(16) is omitted); Gaussian charge distri- 
bution both in nucleus and in proton; exponential distribution for nucleus and 
proton; exponential model for nucleus with proton taken as a point-like object 
(F p =l). With the exception of unrealistic extreme cases, calculation results 
are not very sensitive to the choice of a specific formfactor model. At the same 
time, the difference with the calculations performed for a point-like nucleus 
(dashed curve in the figure) is very big, and reaches the orders of magnitude 
in the probability of large deflections. 

The total average value of squared transverse momentum accumulated 
in the layer can be obtained, if we add the contribution of large deflections 
(q > q b ) to the restricted value given by Eq.(14). For Gaussian charge distri- 
bution, such evaluation can be easily performed analytically: 

W> = 24 h^— + 2Z in ( i + 

where q^, q p , q c , and q a are defined by Eqs.(7)-(9). The average contribution 
of quasielastic processes (the second term in brackets) is small even for light 
elements; for example, in carbon the rms-deviation increases by only 1%. How- 
ever, these interactions have to be taken into account for adequate description 
of the distribution function in the tail region (see Fig.3). 

For three-dimensional calculations, when both angular deviation and lat- 
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eral displacement of muon are of interest, the restricted contribution (of q < q^) 
can also be calculated in a Gaussian approximation, but taking into account 
the well-known correlation between the angle and the displacement [13]. 



4 Atomic electron contribution; mixtures and compounds 

Contribution of incoherent scattering on atomic electrons is often taken 
into account by means of a substitution of Z 2 by Z[Z + 1) in basic relations 
(for differential cross section, some parameters of the theory, etc.). However, 
this substitution is rather inaccurate, since, firstly, the cross section at small 
angles (and low momenta) is determined in this case by the inelastic atomic 
formfactor (with different g-dependence and different characteristic radius of 
the atom). Secondly, due to light electron mass, kinematics of the process (in 
laboratory frame) appreciably changes. That leads to the appearance of the 
upper limit on transverse momentum and on deflection angle, and to an appre- 
ciable probability to loose large fraction of muon energy in a single collision. 
Of course, target size corrections for scattering on electrons are irrelevant. 

Consideration of kinematics of muon scattering on a free electron gives 
the following expression for the transverse momentum q 

q = yj2m e T(l-T/T m ), (18) 

where T and T m are kinetic energy transferred to the electron (initially at 
rest) in the collision and the maximal value of this energy: 

T m = 2m e p 2 /(m 2 e + ml + 2m e E). (19) 

The accumulated squared transverse momentum in the layer X may be eval- 
uated similar to Eq.(13), combining the contribution of \ow-q region (where 
the inelastic atomic formfactor is important) and of high-g domain. In the lat- 
ter case, (Q 2 ) is easily calculated by means of the convolution of the squared 
value of q ( Eq.(18)) and the differential cross section of knock-on electron 
production. In this way, we obtain the following relation: 



(Q 2 ) = 8nN Av jX(r e m e ) 2 

Within the Thomas- Fermi model, K'—1429; an accurate value for hydrogen is 
K'=AA6. Small corrections which correspond to spin terms in the cross section 
of muon-electron scattering are omitted for simplicity. 

Rare catastrophic collisions with large energy transfers are also included in 



In (k'Z~ 2 '\\ 2 ^ 



m. 



(20) 
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(Q 2 ) calculations (Eq.(20)) on the average. For simulation of muon transport 
in thick layers, such collisions are usually treated separately (with explicit 
simulation of kinematic variables). The average accumulated Q 2 in collisions 
with energy transfers less than a certain cut T c (T c < T m ) is given by 



(Q 2 )rstr = 8nN Av jX(r e m e ) 2 

Multiple scattering on the electrons represents only a correction to the scat- 
tering on nuclei, and maximal transverse momenta and deflection angles are 
small. Therefore it seems justified to introduce this correction on the basis of 
a Gaussian approximation (for example, as an addition of contribution given 
by Eq.(21) to the restricted value of accumulated Q 2 on nuclei defined by 
Eq.(13)). Such procedure allows to avoid double counting of the collisions 
with energy transfers T > T c , which will appear if multiple scattering process 
and knock-on electron production are simulated incoherently. 

As a general prescription, for calculation of multiple scattering in mix- 
tures or compounds, the differential cross sections Eq.(4) for elements have 
to be averaged with weights Wi corresponding to relative mass abundance of 
individual entries. 

Within the framework of a mixed Monte Carlo approach described above, 
this leads to: (i) summation of restricted accumulated Q 2 -values with weights 
Wf, (ii) description of the cross section in the intermediate region (q a <C q <C 
<7at) with the average parameter (Z 2 /A) (and, hence, of the integral cross sec- 
tion for q > qb with the same factor); (iii) averaging of the acceptance function 
Eq.(16) for sampling of random momenta q q > g&: 

P(<k) = f {Fm + - / 1 • (22) 

Here, -F/vife) denotes the elastic nuclear formfactor for the i-th element of 
the mixture. For hydrogen, the expression in braces reduces to the first term 
(see also comments after Eqs.(4) and (16)). Similarly, the contribution of the 
scattering on atomic electrons is determined by the weighted sum of partial 
contributions defined by Eqs.(20) or (21), for total or restricted (with T < T c ) 
accumulated Q 2 values, respectively. 



5 Moliere's equations 

In this section, we give the main equations of the Moliere's multiple scat- 
tering theory where the nucleus is treated as a point-like charge. 




.(211 
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Let f(0,t) be spatial-angle distribution function, i.e. for small-angle ap- 
proximation f(9, t) it is the number of charged particles in the angular interval 
of polar angle -j- 9 + d$ after traversing a thickness of t. Moliere solved the 
transport equation for determination of f(0,t) and obtained the general ex- 
pression for angular distribution function as a following: 

oo 

f(0,t) = -JtdtJ o me^\ (23) 
o 



where Jo(£#) is the Bessel function of order 0. For a homogeneous scatterer 
without energy loss, we have 

oo 

n(0 = 2vrf J xd X [J (eO - l]W(x), (24) 



where 

W( X ) = Na( X ). (25) 

Here N = NavP/A (p is density of the scatterer in g/cm 3 ) is the number of 
scattering atoms per unit volume, with differential cross-section cr(x) each. 

The scattering of relativistic particle by an atom is determined by a mod- 
ified Rutherford law 

°(X) = <WxMx)^(x), (26) 



where cr Ru (x) is Rutherford cross-section for scattering by point-like charge, 
*&(x) is a function which takes into account the screening of the nuclear 
Coulomb field by atomic electrons, and ^w(x) is the nuclear formfactor Eq.(16). 
In Moliere's theory, the following approximations are used: JF^(x)=1, 

Z 2 e 4 

^Ru(x) = 4—, (27) 
VX 



and the screening function is 

nx)=xV(x 2 + X 2 J, (28) 

where the screening angle is 

X 2 a = xl(l-l3 + 3.76(Za) 2 ), (29) 
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1.13 



a = 1/137, and xo ^ ——Z 1/3 m e /p. 



137 



In his original paper, Moliere used the variables 
'd = 6/xc^B and r] = t i Xc\ r B, 



(30) 



where the characteristic angle Xc is defined as 
AnNtZ 2 e A 



(31) 



Here t is measured in cm, and B equals to the solution of the transcendental 
equation 



B-\nB = l-2C + \n( X 2 J X l) 



(32) 



(C = 0.5772 is Eiler constant). Then the expression for the angular distribu- 
tion function can be written as 



f{e,t)0M = f M {#,B)-&&&, 



(33) 



where 



2 



Moliere's expansion method is to consider the term [rf ln(r] 2 /A)] /AB as a 
small parameter; then the exponent can be expanded to second order terms 



/m(*, b) = /j? w + fS> w/b + f^m/B 2 + ..., 



(34) 



where 



ri 2 , ri 2 
— In — 

4 4 



(35) 



The first two functions have simple analytical forms: 

/l?(tf) = 2e~* 2 , /j?(tf) = 2e^ 2 (tf 2 - l) [e^ 2 ) - lntf 2 ] - 2 (l - 2< 
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where E 1 ($ 2 ) is the exponential integral. The Moliere's theory is valid for 
B > 4.5 and X 2 C B < 1 [2], [3]. 



6 Modified Moliere theory 



The effect of nuclear size can be evaluated using a suitable formfactor 
3~n(x) hr Eq.(26). The general method to take into account nuclear effects 
for arbitrary charge distribution in nucleus has been suggested by Cooper and 
Rainwater [4] . The Gaussian approximation of nucleus and protons formfactors 
was regarded in paper [5]. In this case, assuming that q p ^> q^ one can re-write 
the Eq.(16) as 



J^( x ) = (i _ Z- 1 ) e-* 2 / fl A +Z ~ 1 e 



xV4 



(37) 



where (see Eq.(7) and Eq.(16)) 

a A = Ml/E 2 and a\ = 3q 2 /E 2 . 

Using Eqs.(24)-(28) and Eq.(37), the function f2(£) can be written as 

X&X 



m = 2 X 



{ (x 2 + xl) 



[MZx) ~ 1] [(! - z ~ l ) e " xV4 +Z- l e-* 2 ' a l] .(38) 



This integral can be evaluated analytically, and we have 

m = - J + \ [(l - Z- 1 ) D A (r), r 2 ) + Z-'D p (r), r 2 ) 



(39) 



where 



D (p, r 2 ) 



r]_ 
4 



In 



E 1 



and 



r 2 = a 2 / X 2 c B. 



Then angular distribution function r) is 

oo 

/ M (0,r) = / lydiyJo^^e^expj- [(l - Z" 1 ) D (??, r 2 ) + Z" 1 !) (77, r 2 )] } .(41) 
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The exponent in Eq.(41) may be expanded in a series, and we obtain 

Im{*, r) = £ j-fV(^ r) » /£>(*) + -ffitf, r) + • • • , (42) 

n=0 B B 



where 



— f r]dr]J {r],'i9)e 

TV. J 




V/4 



(l - Z- 1 ) D (p, r 2 ) + Z- l D fa r 2 )] n .(43) 



At n = and 1, the integrals Eq.(43) equal to Q 

fflp) = 2e"^ , /«(tf,r) = 2e-" 2 [ 2 + (l - Z^) G (^r 2 A ) + Z^G (tf,r 2 



where 



G(V 2 ) = l(c+l + lnr 2 )+/ t 2 



— 1 



,Q2 

E%{--) -In?? 5 



K 2 e^ 2 / K2 ,(45) 



and k 2 = 1 + 1/r 2 (in Ref. [5] the coefficient k 2 is missing in the last term). 
The integral angular distribution P($,t) is 

P (#, t) = £ — / f$ « t)& d& « P (0) (0) + P (1) (0, t) + • • • , (46) 



where 



1 r 



P {0) m = e"" 2 , P (1) (^,r) = ^ [(l - Z- 1 ) Q (tf,r 2 ) + Z~*Q (tf,r 2 )] ,(47) 



and 



2\ „-1? 2 



Q(tf,r 2 )=e 



V + C + l + lnr 2 



+ (k 2 - 1) [R, (#, r 2 ) + R 2 (tf, r 2 )] .(48 



The functions i?i and i?2 are 



Ha (*V 2 ) 



Pi 



?? 2 



1 + r 2 



£i(-tf 2 ) 



1 Note that integrals with E\ (jf'T 2 /4) can be evaluated using derivation over the 
parameter r 2 . 
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R2 r 2 ) 



d 2 



K 2 -\ 



- 1 



Bite)-** 



The function f M ($,r) -> / M (tf) at (R 2 N ) and (i? 2 ) o. On the other 
hand, r 2 — > oo at t — > also. So, the angular distribution of relativistic 
particles should be similar to Moliere distribution after passing through the 
small thicknesses of matter. However, even for thin layers the distribution 
function appreciably deviates from the Moliere law at large deflection angles 
(9 > qw/E), and tends to the cross section behaviour Eq.(26) with formfactors. 
The expansions to the first order terms of Eq.(42) and Eq.(46) are valid up 
to r\ > 0.2. For thick layer of matter (t\ and r 2 <C 1), the assymptotic of 
Eq.(42) is Gaussian [5]: 



OO 

f M (0, T ) ~ / J (rr&)e-^ridri = -e"^, 



(49) 



where 



2-C-lnrl 



T 21 

n 



(50) 



If the mixture of n different scatterers is present, the characteristic angle 
Xc is given by 



Aire 4 

p2 



(51) 



i=i 



i=i 



and the integral angular distribution function t) is defined by 



p(#,t) = Y,^z 2 p^,t) / ^WiZ?, 

1=1 / 1=1 



(52) 



where the subscript % denotes the atomic species and weight w,i is proportional 
to the number of atoms of the i-th element in the mixture. 

The integral angular distributions on nuclei as a function of the transverse 
momentum Q calculated by the mixed Monte Carlo tecnhique and analytically 
Eq.(46) are shown in Fig. 4 (iron targets of 1 cm and 100 cm) and in Fig.5 (wa- 
ter thickness of 1000 cm). The results are in good agreement. Moliere angular 
distributions are shown also for comparison. A minor difference between ana- 
lytical and Monte Carlo calculations in the tails of the curves is conditioned by 
a different choice of the proton formfactor (Gaussian charge distribution with 
R p = 0.80 fm in analytical evaluation and exponential one with R p = 0.85 fm 
in Monte Carlo simulation. 
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7 Conclusions 



A high-energy muon is capable to cross very thick layers of material and its 
de Broglie wave is comparable with nuclear dimensions. Hence, the probability 
of scatters with impact parameter of the order of nuclear radius ( scatters with 
large angles) is not small. For such interactions, the law of single scattering 
has to be modifed significantly to take into account the finite size of charge 
distribution in the nuclei. Nuclear size effects can be desribed by multiplying 
the Rutherford cross section by a nuclear formfactor JF^, which goes to 1 for 
large impact parameters and vanished at very large transverse momenta. 

We have presented two methods for calculating the multiple scattering 
distribution for muons traversing a thick scatterer. The mixed Monte Carlo 
technique and modified Moliere's theory of MCS take into account the nuclear 
size effect. Parametrisation of A-dependence of rms-radii for light and medium 
nuclei is given. The results of two methods are in good agreement, and show 
that the muon angular distribution is similar to Moliere distribution after 
passing through small thickness of material and moderate deflection angles. 
However at large thickness of scatterer, the angular distribution is drastically 
changed by the influence of the nuclear formfactor. So, in this case the Moliere 
theory heavily overestimates the probability of large angle deflections. Con- 
tribution of incoherent scattering on atomic electrons is also correctly taken 
into account. 

A necessity of taking into account nucleus size effect as well as incoherent 
scattering on atomic electrons when describing the muon scattering at large 
angles in thick matter layers is shown. Neither Gaussian approximation of 
multiple scattering distribution nor Moliere theory can be used for accurate 
description of high-energy muon transport in a wide range of the thicknesses 
of material. 
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Figure Captions 



Fig.l: Differential distribution in accumulated transverse momentum for muons 
after passing iron of 1 cm. Dashed curves, open circles: point-like nucleus. Solid 
curves, dark points: calculations taking into account the nuclear formfactor. 
Curves are calculated with g&=l MeV (10 8 events), the points correspond to 
q b =0.1 MeV (10 6 events). 



Fig. 2: The same as in Fig.l but for iron target of 100 cm. In this case, the 
values of q b in simulation were 10 MeV and 1 MeV, respectively. 



Fig. 3: Integral distributions in the value of accumulated transverse momentum 
for muons after passing iron layer of 100 cm, calculated with different models 
of the nucleus (see the text). 



Fig. 4: Integral distributions in the value of accumulated transverse momentum 
for muons after passing iron layers of 1 cm and 100 cm calculated by the mixed 
Monte Carlo tecnique (open circles for point-like nucleus and dark points for 
finite nuclei) and analytically (solid curves are result of Moliere theory and 
dashed curves are calculations taking into account the nuclear size effect). 



Fig. 5: Integral distributions in the value of accumulated transverse momentum 
for muons after passing water layer of 1000 cm calculated by the mixed Monte 
Carlo tecnique (dark points for point-like nucleus, open circles for finite nuclei 
and stars is a result of calculations, taking into account the nuclear formfactor 
and electron contribution) and analytically (solid curves are result of Moliere 
theory and dashed curves are calculations taking into account the nuclear size 
effect). 
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